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ABSTRACT 



Preliminary Turbine Design for Thermal Cycle Feasibility 

David A. Potter 

Submitted to the Department of Ocean Engineering and the 
Department of Mechanical Engineering on 10 May 1974 in partial 
fulfillment of the requirements for the degree of Master of 
Science in Naval Architecture and Marine Engineering and 
Master of Science in Mechanical Engineering. 

Mathematical models of an axial flow turbine and a radial 
inflow turbine were developed. Utilizing these models, a com- 
puter program was written such that given a thermal cycle and 
a minimum of input variables, the necessary turbine design cal- 
culations are performed. An extensive output is provided such 
that the designer may utilize this output in order to evaluate 
the practical feasibility of the given thermal cycle. 

The Zener sea cycle is used as a test case with axial and 
radial designs examined for feasibility. 

The computer program is written in the Fortran IV language 
specifically to be run on the Interdata 70 computer. 

Thesis Supervisor: A. Douglas Carmichael 

Title: Professor of Power Engineering 

Thesis Reader: Warren M. Rohsenow 

Title: Professor of Mechanical Engineering 
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NOTATION 



Blade cross sectional hub area 

Flow area 

Annulus 

Blade chord length 

Specific heat at constant 
pressure 

Specific heat at constant 
volume 

Hydraulic diameter 
Gravitational constant 
Enthalpy 

Moment of Inertia in 
tangential direction 

Moment of Inertia in 
Axial direction 

Conversion factor 

Windage loss 

Stator loss 

Rotor loss 

Clearance loss 

Mass flow rate 

Mach Number 

Bending moment in tangential 
direction 

Bending moment in axial 
direction 



UNITS 

ft 2 

ft 2 

ft 2 

ft 2 

BTU/°F 



BTU/°F 



ft 

32.174 ft lbm/lbf sec" 
BTU/lbm 
ft 4 



ft 



778 ft-lbf/BTU 
BTU/lbm 
BTU/lbm 
BTU/lbm 
BTU/lbm 
lbm/ sec 
Non-Dim 
f t-lbf 



ft-lbf 



0 

p 

e 

R 

a 

r 

R 

e 

RPM 

S 

s 

T 

t 

U 

•V 

V 

x 

ex 

/3 
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lo 

I, 



Number of axial stages 
Specific speed 
Throat opening 
Pressure 

Volumetric flow rate 
Universal gas constant 
Reaction 
Radius 

Reynold's Number 

Shaft rotative speed 

Spacing 

Entropy 

Temperature 

Blade height 

Rotor tangential speed 

True fluid speed 

Relative fluid speed 

Axial distance 

True fluid flow angle from 
axial 

Relative fluid flow angle 
from axial 

Ratio of C /C 
P v 

Primary loss factor 

Primary loss factor corrected 
for aspect ratio 



Non-Dim 
Non-Dim 
feet 
lbf/in 2 
ft"V sec 
13TU/lbrn °R 
Non-Dim 
feet 
Non-Dim 
Rev/ min 
feet 

BTU/lbm °R 
Deg Rankine 
feet 
ft/ sec 
ft/ sec 
ft/ sec 
feet 
Peg 

Deg 

Non-Dim 

Non-Dim 



Non— Dim 
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? Primary loss factor corrected 

f 3- for AH & R 

e 

^ Slip factor 

yj Efficiency 

© Tangential angle 

^ Stage loading coefficient 

JU Static viscosity 

^ Kinematic viscosity 

A Differential function 

<T* Solidity 

<5^ Centrifugal stress 

c 

VT Constant 

jO Densi ty 

<j) Plow coefficient 

ty Work coefficient 

Angular rotative speed 



SUBSCRIPTS 



Non-Dim 

Non-Dim 

Non-Dim 

Deg 

Non-Dim 
lbm/hr ft 
ft 2 /hr 
Non-Dim 
Non-Dim 
lbf / in^ 
3.1415927 
lbm/ft 3 
Non-Dim 
Non-Dim 
Radians/ sec 



0 Stagnation conditions 

01 Machine inlet stagnation condition 

02 Machine outlet stagnation condition 

o2s Isentropic outlet stagnation condition 

i Individual flow path 

t Total flow paths 

B Tangential direction 



X 



Axial or meridional direction 



m 


Mean 


h 


Hub 


t 


Tip 


f 


Fluid 


tt 


Total to total 


ts 


Total to static 


mat. 


Material 


max 


Maximum 



Axial 



1 


Rotor inlet 


2 


Rotor outlet 



Radial 



1 


Stator inlet 


2 


Stator outlet 


3 


Rotor inlet 



4 



Rotor outlet 
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INTRODUCTION 

The thermal cycle is the basic tool in use today to con- 
vert the energy contained in a source or fuel into useful 
power* Cycles using geothermal energy, cycles powered by 
solar energy, and cycles powered by the thermocline existing 
in the oceans, are just a few of the proposed cycles being 
studied as possible sources of power for the future. In the 
search for new and plentiful energy sources in the years to 
come, new and different thermal cycles will be proposed. 

The evaluation of a given thermal cycle as to its prac- 
tical feasibility involves a detailed study in which not only 
the energy source, but the hardware needed to convert that 
energy into useful power, are investigated. In such a study, 
the engineer must concern himself with three main areas: con- 
verting the energy source to heat, transferring that heat 
energy to a working medium, and then utilizing the increased 
energy of the working medium to produce useful work or power. 
It is with this last facet of the problem that this thesis is 



concerned . 
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BACKGROUND 

1,1 Thermal Cycles 

Thermal cycles come in all sizes and shapes. They are 
used not only to convert heat energy to useful work, but to 
transfer heat energy from one location to another, as in a 
refrigeration unit. The difference between the two cycles, 
other than hardware, lies in the fact that in a work-orienta- 
ted thermal cycle energy is supplied in the form of heat or 
fuel and work is the end result; while in a heat transfer 
or refrigeration cycle, work is supplied as an input and the 
end result is the transfer of heat energy. While both types 
of cycles operate under the same principles, it is the former 
type which is of great concern in the search for means to 
supply society’s growing energy demands, and which will be 
the only one considered from hereon. 

Thermal cycles of the heat to work type can be catego- 
rized as either open or closed cycles. In the open cycle, 
the energy of the source (fuel) is transferred to a working 
medium in the form of heat, whereupon this increased energy 
of the working medium is utilized to produce useful work. 

At the conclusion of the open cycle, the working medium is 
exhausted from the machine. In a closed cycle the process 
is essentially the same; however, the working medium is re- 
tained in the cycle after producing work and is reintroduced 
to the source where it again increases its energy state with 
the addition of heat. There are, of course, many variations 
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and combinations of the above cycles using many different 
working mediums at many different temperatures and pressures. 
1,2 Turbomachinery 

One problem common to all thermal cycles, open or closed, 
is that of finding an efficient means to convert the energy 
of the working medium into useful work. Without considering 
unproven schemes, such as magnetohydrodynamics, there are 
two basic classes of machines in use today with which this 
can be accomplished: machines which extract the energy from 

the working medium as a change in its potential energy, and 
those machines which extract the energy as a change in the 
kinetic energy of the working medium. Without discussing 
the relative merits of either type of machine, it is suffi- 
cient to state that the second class of machine, the turbo- 
machine, is a viable alternative and one which should be 
considered in any thermal cycle feasibility study. 

Turbomachines, or more specifically, turbines used in 
thermal cycles today, can be divided into two categories: 
those in which the working fluid flows through the machine 
in a nearly constant axial direction parallel to the center- 
line of the machine; and those in which the working fluid 
flows in a direction radial to the centerline of the machine, 
either inward or outward. The choice as to which type of 
machine is best for a certain set of conditions is a diffi- 
cult one, depending on a large number of variables external 
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to the cycle itself, such as weight, cost, and operating 
load. This will be discussed more fully in following 
sections. 

1.3 The Design Sequence 

In determining the technical feasibility and practical 
value of a given cycle, at the preliminary design stage, the 
engineer must follow an iterative process. One logical 
sequence would be to: 

1) Assume a working medium and a set of states for 
the medium (most commonly the energy extraction 
process is done in a gaseous state). 

2) Assume certain hardware restraints and a flow 
rate for the working fluid. 

3) Design the required hardware (if not available 
commercially) . 

4; Evaluate the necessary hardware and if not 
satisfactory return to step 1. 

This process is repeated several times until the hardware 
conforms to the engineer's definition of feasibility, or until 
some prechosen parameters are optimized. The problem here 
lies in the fact that the calculations necessary to design the 
hardware, although relatively straightforward, are numerous 
and time-consuming. For a given set of cycle conditions, 
there are often many options open to the designer, all of which 
will result in different .hardware and all of which will meet 
the design criteria. The choices between axial or radial 
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turbines, and between turbines or reciprocating machines, 
are good examples of this. There is a need, therefore, for 
the designer or engineer to quickly translate his cycle re- 
quirements into specific hardware sizes and characteristics 
in order to logically choose between the options available, 

1.4 Optimization 

Feasibility and optimization are two words whose 
definitions can be different for every different cycle con- 
sidered, In a given situation it may be necessary to max- 
imize efficiency, while in another situation, the design 
may be constrained such that weight or volume or cost must 
be minimized. As for turbines, in some applications an 
axial flow unit might be chosen for its high efficiency, 
while in others a radial machine with its lower fabrication 
cost might be best. Feasibility also has different meanings 
in most different contexts. In one instance it might mean 
technically possible, while in another it might mean finan- 
cially profitable. In order to answer these questions, the 
engineer or designer must first define feasibility and then 
decide upon which, if any, parameters are to be optimized 
within the feasibility constraints. Having done this, he 
may then evaluate the hardware designs and make a logical 
choice between competing designs. 

i 

1 . 5 Purpose 

Within this framework then, the purpose of this thesis 
is to bring together the various concepts of turbine design; 
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to develop a mathematical model; and to computerize that 
model in such a way that, given a minimal amount of in- 
formation about the cycle and the machine constraints, a 
set of turbine hardware parameters will be provided that 
meet the requirements of the preliminary design stage. 

These may then be used in the design sequence mentioned 
previously to determine feasibility. No attempt will be 
made to optimize within the computer program, but rather 
to provide information from which a suitable optimization 
criteria can be evaluated. For this reason no attempt will 
be made to choose between an axial or a radial turbine. 

Both will be developed simultaneously and information will 
be provided on both types. The program will be capable of 
handling multiple inputs for the same cycle in one computer 
run . 

It is hoped that by utilizing this program the design 
sequence may be shortened considerably. This will be true 
especially if use is made of other computer programs to 
calculate heat exchanger parameters. In doing so, the 
designer will be able to study more completely the possibi- 
lities of a given cycle at the preliminary design stage, 
and thus develop a more suitable final design. 
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DEVELOPMENT OF THE MATHEMATICAL MODEL 

2.1 General 

A design philosophy is a necessary ingredient to any 
worthwhile design. This philosophy may not always be ex- 
plicitly stated; however, it is always present in a good 
design. It is a means whereby the designer is able to make 
logical and, of most importance, consistent decisions during 
the course of the design sequence. This is especially true 
at the preliminary design stage where decisions are made 
which constrain the course of the design throughout. 

The philosophy followed in compounding the mathematical 
models described in the following sections and in the actual 
programming itself can be stated as accurate and simple. 

Both of these qualities apply to any work done at a prelimi- 
nary design stage. Accuracy is required because the designs 
have now passed beyond the back of the envelope stage; and, 
although detailed design is not required yet,. the results 
must be accurate enough so that detailed design can begin by 
using them as a basis. Simplicity is a must because at the 
preliminary stage, neither the time nor the money are usually 
available to carry out the extensive calculations of a de- 
tailed design. That is not the purpose of the preliminary 
design. One last idea is also followed in the handling of 
certain relationships; that in the absence of accurate and 
simple relationships to cover a certain situation, traditional 
and reasonable approximations are assumed. With the lack of 
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specific data in certain areas this borders on an educated 
guess, but then again this is what the designer is paid for. 

As to the development of the models themselves, they 
can be divided into three main areas: representation of the 
working medium itself, representation of the axial turbine, 
and representation of the radial turbine, They will be ex- 
plained in that order. 

It is assumed that the machines described by the models 
can be built of conventional materials and that the axial 
blading has no twist and is set in a purely radial direction. 
It is further assumed that internally cooled blades or rotors 
are not utilized. 

All computations are in English units with the notation 
as explained on pages 8 to 11, Calculations are carried 
through to a precision consistent with the accuracy of the 
Interdata 70 computer (six significant figures). 

2.2 The Working Medium 
2.2.1 Equation of State 

Two assumptions concerning the working medium are basic 
but should be stated for clarity. It is assumed that the 
working medium is a fluid and that it enters the machine in 
a gaseous state; this precludes the design of hydraulic tur- 
bines. The second assumption is that the inlet and exit 
states of the fluid arc* known and that the designer has the 
properties of the fluid at these states available as input. 

In order to accurately represent the fluid within the 
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machine, an equation of state for the fluid had to be picked 
The perfect gas model was chosen over the other equations of 
state for its flexibility, its accuracy, its vide acceptance 
and its compatibility with the standard one dimensional com- 
pressible flow relationships. It is therefore assumed that 
in the region under consideration^ the fluid can be repre- 
sented by the perfect gas model and that the effects of 
velocity can be represented by the standard one dimensional 
compressible flow relationships. 

To make this assumption more realistic, it is necessary 
to calculate the different constants in the region under con 
sideration rather than use traditional or tabulated values. 
One method to achieve this is to use data for an isentropic 
expansion in the region and to calculate the constants from 
this. Since the designer already has calculated the point 
in order to determine his outlet conditions, it was decided 
to use an isentropic expansion across the machine for the 
reference points. Thus from points 1 to 2s, 

C = (h , -h , )/(T 0 -T , ) 
p o2s ol o2s ol 

Using the Gibbs equation and the perfect gas law, the 
following commonly used relationship for an isentropic 
expansion may be derived: 

(r 7 /p , ) = (t ~ /t ) * 

o2s ol o2s ol 



where -C / C 
p v 
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Solving this equation for Ji 1 : 





Thus C = C / ,a 
v p 7 * 

and R = C - C . 

P v 

2.2.2 Flow Rate 

In describing a particular turbomachine, the mass flow 
rate of the machine has a large effect on its design. For a 
particular cycle, the mass flow rate to achieve a given 
output can be easily calculated using the steady flow energy 
equation. It may not be advantageous, however, to achieve 
this output on only one machine, but rather to utilize several 
machines, connected in parallel, to achieve the given output. 
The mass flow rate through a given machine, with no allow- 
ance for leakage, will be: 



2.2.3 Compressible Flow Relations 

The two parameters necessitating the use of the one 
dimensional flow relations are the calculation of Mach number 
and the calculation of flow area. Using these relationships, 
and given the stagnation or total temperature and the 
velocity of a fluid at a point, the Mach number and flow 
area may be easily calculated. From the definition of total 
temperature and Mach number: 



ra 



individual 



total 



/ # units in parallel 



T = T - V 2 / 2g JO 



o o p 
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M = V/(^g o ItT)» 

In its common form: 

m.(T o )VP 0 A f = M(^g o /R)V(l + (y-l/2)M 2 ) (y+l) / 2( ^" 1 ) 

Thus A f = (m i (T o )'VP 0 ) x f(y, g() ,R,M) 

(T ) ^/P Q ) x (flow function) 

where : 

Flow function = (l + l<P-l/2)M^ ) (^-1 )/^j (^g^/K) 2 

2.3 Axial Turbine 
2 . 3.1 General 

The first step in the design sequence involves deciding 
the general configuration of the machine. Is the hub radius, 
the mean radius, the blade tip radius or some combination of 
these to be maintained constant throughout the machine? 

Wilson (l) recommends that the mean radius be maintained 
constant, while Carmichael (2) recommends hub radius be kept 
constant at the preliminary design stage. Constant mean 
radius was chosen for the reason that the effect on turbine 
performance between the two is small, and calculations at a 
later stage will be simplified with a constant mean radius. 
Following this assumption, mean conditions are calculated 
first. As the velocity triangles and other blade parameters 
vary by stage throughout the machine and as the extremes are 
represented by the first and last stages, only the first and 
last stages will be calculated in detail with the other 



stages assumed to lie in between these values. Finally, 
machine parameters which depend upon the previous calcu- 
lation, such as total efficiency and machine length, will 
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be made. It is also assumed at this point that the machine 



is a full admission turbine operated under a constant load; 
both assumptions are consistent with stationary power plant 
operation. 

2.3.2 Required Output 

By a process of elimination, the list of quantities 
which is sufficient to describe the turbine at the prelimi- 
nary design stage was narrowed to the following. Many other 
quantities could have been included; however, all major para- 
meters are included, and those omitted may simply be calcu- 
lated from those given. 

1) First and Last stage Tip Diameters 



2 ) 

3) 

4) 

5) 

6 ) 



II 



It 



If 



II 



II 



Blade Height 



Aspect Radios 
Number of Blades 
Efficiency ( Total to total) 
Velocity Triangles 
Critical Mach number 
Efficiency ( Total to total) 



7 ) Machine 



8 ) 

9) 

10 ) 



II 



II 



fl 



Critical Stresses 



Length 
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In describing the velocity triangles for a given stage, 
three locations for both the stator and rotor entrance tri- 
angles are described: hub, mean and tip. 

2.3.3 Input 

In listing the factors which are generally used to des- 
cribe a turbine stage, it was found that they could be cate- 
gorized in three groups: those which describe the mass flow 
rate; those which describe the work accomplished per pound 
of fluid per stage; and those which describe the performance 
of a given stage. Each facet is necessary for the preliminary 
design . 

The mass flow rate for the machine, and thus each stage, 
(since the stages are assumed to be in series with no bleed- 
offs) can be described by the total mass flow rate for the 
cycle and the number of turbines in parallel, as described 
in section 2.2.2. 

The work per pound of fluid per stage can be simply re- 
presented also. Making the assumption that the work accom- 
plished per stage is equal for all stages (2), the individual 
stage work can be represented by the total work or enthalpy 
drop for the machine and the number of stages. 

The parameters used to describe stage performance are 
not so straightforward. They can be classed as either pri- 
mary or secondary, depending upon the strength of their in- 
fluence upon performance. Primary factors are generally 
associated with the geometry of the velocity triangles, and 
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if the deviation is assumed negligible, the geometry of the 
blade itself. Secondary factors have a lesser effect on 
efficiency and are related to the likes of clearance, aspect 
ratio, and other mainly mechanical factors. In order to main 
tain a simple machine only primary factors will be specified 
as an input, with a logic developed for the determination of 
reasonable secondary factors based on the given primary 
factors . 

As primary factors, Carmichael (2) recommends reaction, 
and inlet flow angle, c< . Craig and Cox (3) prefer to 
use the stage loading parameter, V, also called the work 
coefficient, and the flow coefficient,*^. 7 \ is another widely 
used parameter (4). Comparing the definitions of these 
parameters , 

ft = 1 - * v ei + v e2 >/ 2lJ 

Tan (oCj) = V 0 /V x 

V =iuv e /u 2 
= v x /u 

*• = u /< v ei - v 02> 

it can be seen that there are only four independent variables 
For reasons explained in a later section, one of these can be 
eliminated; thus only three independent variables will des- 
cribe the shape of the velocity triangles. 
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The three variables chosen for use here are~R!,U , and ^ . 
Reaction was chosen for its compatibility with the Euler work 
equation which will be used later to solve for a mean radius. 
Given this mean radius and an KPM, the second variable, U, 
can be specified. The last var iable , oC ^ , was chosen as it 
would have to be calculated anyway during the course of the 
design. 

The input for the axial turbine can be summarized as: 

1) Mass flow rate 

2) Number of units in parallel 

3) Number of stages per unit 

A) Shaft revolutions per minute 

5) Reaction at mean conditions 

6) oC^ at mean conditions 

7) Cycle conditions of section 2.2.1, 

2.3.4 Mean Conditions 

The starting point for the mean conditions is with the 
simultaneous solution of the Euler work equation and the 
reaction equation. It is assumed that the fluid enters the 
stator only in an axial direction, with no tangential compo- 
nent. This is a common assumption (4) and is made for several 
reasons. One argument is that the turbine works by extracting 
energy from the fluid; thus any energy left in the fluid at 
rotor exit (stator inlet) in the form of vorticity (tangential 
velocity) is unavailable for extraction, and thus contributes 
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to a loss. Using this assumption, the Euler equation applied 

to mean conditions across the rotor becomes 

= to r V /g J 
o m 01m' & o 



n 



where Ah = h . -h 0 
o o 1 o2 

The reaction equation, also applied at mean conditions across 
the rotor is simplified to 

71 = 1 - V„, /2wr . 

m 0 1 in m 

The simultaneous solution of these equations for r and 

n m 01m 



with Ci)=7fRPh/30 gives 



m 



= (Ah g J450/rr 2 RPM 2 n (l-1 : l) )* 

oo m 



©lm = r -Tr'RPM (l-7l)/l5. 
m m 

Substituting in the appropriate constants, 

1 

r = 1068. 247 (Ah /n(l-"R. ) ) 2 /liPM 
m o m 



and 



V elm - 223.7331 Uh 0 (l-R)/n)*. 



In order to prevent the solution from expanding at T^ m = I > 

must be limited to less than 0.99. 

The rotor speed can be calculated as 

U. = U„ = U = VrRPM r / 30 . 
lm 2m m m 

Using the tri gnoraetr ic relationships, Figure 1, and the 



c<^ which was given as an input, the rest of the velocity 
triangle for the rotor inlet can be calculated. 
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Figure 1 

V, = v m /Sin(oC ) 

lm 91m' ' lm 

V = V_. /Tan (c< ) 

x 91m' lm' 

/& lm = Arc Tan (Tan (o^ m ) (1—1/ (2— 2T^) ) ) 

W. = V /Cos (A ). 

lm x lm 

Making the assumption that the axial velocity is constant 
throughout the machine as was stated before, the tangential 
velocity at rotor outlet is zero; thus the velocity triangle 
for the stator inlet (rotor outlet) can be calculated: 

= ArcTan(Tan(cc )/2(l-TL)) 

' 2in lm in 



W 2r = V 81m / ^ os ^2m )Tan 
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2m 



V 



02m 



0 



V 



xl 



V 



x2 



V 



x 



V 



2m 



V 



x 



The total gas path deflection for the rotor is the sum 
of the inlet and outlet angles (/^ + ; while for the 

stator it is just , as the inlet angle is zero* 

2*3.5 Stage Conditions 

In solving for stage conditions it can be assumed that 
the extreme values are represented by the first and last 
stages, with all other stages falling in between. If the 
variance is assumed linear, the dimensions and velocities 
can be calculated for any stage with reasonable accuracy. 

This assumption also allows the calculation of machine 
efficiency as the mean of first and last stage efficiencies. 
Based on these assumptions, the following calculations can 
be made for the first and last stages. 

Using the one dimensional compressible flow relationships 
described in section 2.2.3, the flow area can be calculated for 
a given velocity and stagnation temperature. The flow area 
and annulus area are related as shown in Figure 2. 



Using the definition of annulus area and the mean radius, 
the blade height may be found. 



Area 



ann 



Area^T /Cos(oc) 
flow' 




Figure 2 



Area ann =Vri ' tr t 2 - r h 2) 

= ^^t- r h )(r t +r h> 



=Yft2r 



ID 



Blade height = t = Area /2 r It 
6 ann m 



whe 



re r = (r +r )/2 
m t h 



The tip and hub radii are found by respectively adding or 
subtracting one half of the blade height to the moan radius. 

The next step involves making the assumption that a 
streamline of fluid has no appreciable curvature in the radial 
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direction as it moves through the machine. This is not com- 
pletely correct, but for small wall curvatures it can be con- 
sidered accurate (l). This assumption means that a segment 
of fluid, as shown in Figure 3, must be in radial equilibrium. 



RADIAL EQUILIBRIUM FORCES 

p + dp dr 
dr 




Figure 3 

Summuation of the forces acting on this segment gives the 
equation for this condition of radial equilibrium. 

i- 4* = V e 2/g A 

/f dr 

Using the Gibbs equation and the definition of total 
enthalpy for an adiabatic condition, the radial equilibrium 



32 



equation can be rewritten as 



dh - Tds = 1 
o 



dV + 1 
x 



<iU 2 V e 2 ). 



dr dr 2g dr 2g r dr 

& o & o 

iMaking the assumption that the enthalpy drop is uniformly 

distributed over the blade radially, and that entropy is not a 

function of radius, the result is 

o O O 

dV = -1 d v ) . 

x 2 9 

dr r^ dr 

At this point one of several distributions for V or 
may be chosen; however, having made the assumption previously 
that axial velocity is constant throughout the machine, the 
above reduces to 



r Vq = constant. 

This constant may be calculated from the mean conditions 
and, using the above relationship, the tangential velocities 
at hub and tip can be calculated. The blade speed is easily 
calculated at hub and tip as shown previously; and, using 
trignometric relationships, the complete velocity triangles 
for rotor and stator are easily calculated in the manner 
illustrated in section 2.3.3. 

The choice of aspect ratio (blade height/blade chord) 
must now be made. As this is a secondary variable, section 
2.3.3, a means must be derived for a logical decision. A 
study of data presented by Forrester (5) shows that the 
secondary losses attributed to aspect ratio effects decrease 
with an increasing aspect ratio; thus a large aspect ratio 
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would be advantageous from an efficiency standpoint. On the 
other hand, Horlock (6) states that the bending stresses on 
the blade may place a limit, on the blade chord length and 
thus limit the feasible aspect ratio. Any aspect ratio 
chosen should therefore be as large as possible without ex- 
ceeding a predetermined stress limit. Carmichael recommends 
a limit of one ton per square inch which was used in the 
following derivation. 

In order to simplify the calculations necessary, the 
fluid was modeled as incompressible and the blade shape was 
simplified so that the resulting criteria could be applied to 
the general case. The maximum stress was assumed to occur in 
the rotor at the root of the blade and located on the 



Blade Root Stresses 





Figure 4 
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upstream side. By considering the change in tangential momen- 
tum across the blade, the bending moment in the tangential 
direction (Figure 4) can be. shown to be 

M e =/VfV r(r " r h )( V dr 

r h z 

where z = number of blades per stage* 



Integrating this 


equation, and substituting the Euler 


work equation and the 


definition of blade speed gives 


M - 

e 


mesh J(r, -r, )^60 S-, 
o t h 1 

(r t 2 -r h 2 ) (2rr) 2 RPM r h n 



Simplifying and substituting in the constants yields 



11 

CD 

s 


591.209 nifth (r.-r, ) S 1 /RPM r r n. 

o t h 1 m h 


Since Tan (/? ) = 
r m 


Forces^ / Forces , . , . 

0 direction x direction 



the bending moment in the axial direction can be described by 



M = 

X 


M„ Tan ( ,0 ) . 
0 m 


For a particular 


design let Kl be a constant representing 


Kl = 


mall (r. -r, ) 591.209/RPM r, r n 
o t h h m 


then M 0 = 


Kl S 1 


• M = 

X 


Kl S, Tan (A ) . 
1 'm 



For the blade shape shown in Figure 4, the moments of 
inertia about the respective neutral axes can be shown to be 



M 

CDI 

II 


c S 1 3 /69.429 


i- = 

X 


Sj c 3 /72 . 
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For stress at the point A, 

0/1- = 30 . 86/ c S 1 2 

x/lg = 36/S lC 2 . 

The bending stress is then 

cr b = Kl(Tan(^)36/c 2 + 30.86/c S ;L ) . 

Using Zweifel’s criteria to relate S and c (4), 

0.8 = 2 Cos 2 (^ 2h )(Tan(^ h ) + Tan(^ h )) S 1 /c 

substituting in the bending stress equation and solving for 
chord, using the one ton per square inch stress limit men- 
tioned previously, 

c = (K2 f)2 

where K2 = K1 32 . 174/ (2000 )( 144 ) 
f = 36 Tan(^ lm ) + 77.15 Cos 2 (/? 2h ) x 

(Tan(^ lh ) + Tan(/? 2h )). 

This approximation for chord is based on a stress limit 
caused by gas bending forces; however, to account for the 
fact that other design conditions may impose a lower limit 
than this criteria, a maximum limit to aspect ratio of 10 
will be imposed. 

Zweifel’s criteria can now be used to determine spacing 
for the rotor from this calculated chord. Assuming the same 
chord for the stator, spacing can be calculated there also. 
As shown in Figure 5, the throat opening can be figured for 
both the rotor and stator. 
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Throat Opening 



Axis of machine 




.FifTure 5 

°i = s i 

°2 = S 2 t: ° s(/< 2in ) 

Knowing the spacing, the number of blades for rotor or 
stator is simply 

Z = 2'rt'r,/ S. 
h 

•> 

Using the hydraulic diameter of the throat opening, 
again for both rotor and stator, 

D = 2 0t/(0 + t) 

the Reynolds number for both the rotor and stator can be 
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calculated as 



“ e 2 = U 2 



ltc l = U 1 V ml^- 

2.3.6 Efficiency 

The determination of efficiency for axial turbines is a 
much-debated subject, Craig & Cox (3), Smith (5), and Ainley 
and Matheson (5), all correlate their data using the work 
and flow coefficients mentioned previously; while Soderberg 
(4) correlates using gas deflection angle. In either case 
the primary loss is related to blade profile loss. Secondary 
losses are attributed to many factors, some of them being: 



1) 


Clearance lo 


ss 


•2) 


Reynolds number effects 


3 ) 


Aspect ratio 




4) 


Guide loss 




5) 


Runner loss 




6) 


Guide gland 


loss 


7) 


Lacing wire 


loss 


8) 


Balance hole 


loss 


9) 


Wetness loss 




10) 


Disc windage 


loss 


11) 


Partial admi 


ssion loss 


12) 


Annulus wall 


loss 


13) 


Cavity loss. 





While including all of these factors would give a most 
accurate efficiency prediction, the accuracy required in a 
preliminary design doesn't require the inclusion of them all. 

Soderberg’ s correlation was chosen in that it provides 
accuracy consistent with the intent of the program; yet it 
is simple enough to be easily programmed on a computer. This 
correlation takes into account profile loss, clearance loss, 
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aspect ratio, and a Reynolds number effect. These losses are 
equated to a fraction of the kinetic energy available to the 
turbine and a total to total efficiency calculated. 

The method consists of finding a primary loss factor, 
based on gas deflection angle for both stator and rotor. 

This loss factor is modified for aspect ratio and Reynolds 
number effects, and then an efficiency calculated. The 
efficiency is then modified by a factor which takes into 
account the clearance losses. 



The above are calculated separately for both stator and rotor. 
The efficiency is then: 



The primary loss coefficient is one which can be pro- 
grammed easily in a curve-fitting subprogram. The tip clear- 
ance area was calculated based on a tip clearance of one per- 
cent of the tip radius. This gave a constant factor of 0.99. 
Assuming a linear distribution of efficiencies between the 
first and last stages means that the machine efficiency is 
the average efficiency. 



? = f (Gas Deflection) 

* o 

? 1 = ((l+?J (.975+0.075 c/t))-l 





01m 



>^ tt (Stage) = (A 



rea Annulus - Tip Clearance Area ) . 
Area Annulus 
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^ Machine ~^J?irst Stage +> ^Last Stage * 
2.3.7 Miscellaneous Machine Parameters 

In order to more easily picture the machine size, all 
pertinent radii are converted to diameters* 

A machine unit length is calculated using the average 
stage chord and the number of stages. 

The maximum blade centrifugal stress occurs in the last 
stage at the rotor blade root. This can be written as 
follows : 



<r 

c 



= nat . 

g A u 
& o h 




2 

co r 



A u dr . 
h 



If the blade is assumed to have no taper, A^ is constant. 

Substituting in the constants and assuming the density of 

steel is approximately 0.29 lbm/in gives 

<y = 0.001186 RPM 2 r t. 
c m 

2.3.8 Turbine Limitations 

Certain limitations must be placed on the mathematical 
model previously described, so that it won’t produce turbine 
parameters for a machine that is either physically impossible 
or not consistent with the assumptions of the model. These 
limitations will be in the form of numerical checks in the 
program and will cause the calculations to be terminated if 
they are not satisfied. 

The first check is for excessive gas deflection. While 
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this in itself is not an impossible situation, only leading to 
a lower efficiency, the data available for the Soderberg 
correlation only accommodates a gas deflection of up to 140°. 
Since the efficiency decreases rapidly at angles greater than 
this, no harm is done in limiting gas deflection to 140°. 

Reaction is another parameter which is limited. The 
upper bound has been discussed earlier in section 2.3.3. A 
check must be made, however, to insure that at no point does 
reaction become negative. Physically, a negative reaction 
would mean that within a given stage, at some point, a pressure 
rise was occuring to be compensated for by a larger pressure 
drop at another point. The overall condition, of course, 
would be a pressure drop. In order to avoid this condition 
which brings about higher losses, a check should be made to 
insure non-negative reaction. Negative reaction can only 
occur with large Vg and since it is assumed that r Vg = 
constant, the check should be made at the smallest radius, 
the last stage hub. 

In solving for a mean radius, in section 2.3.4, no con- 
tinuity check was made. without a check for continuity it 
would be possible, using the model described, to compute a 
negative machine diameter. In order to avoid this and to 
allow for bearings, shafting, and blade attachment, a check 
will be made to insure that the hub radius calculated is 
greater than 50 percent of the tip radius. 
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The one dimensional compressible flow equations used 
throughout also pose another limitation. The fluid flow 
throughout the machine is assumed to be subsonic, and 
the design has proceeded such that no consideration has been 
given to the shock effects caused by supersonic flow. To 
insure subsonic flow, a check will be made on the two largest 
velocities to insure a Mach number less than one. The 
velocities are the last stage hub velocity at rotor inlet and 
the last stage tip relative velocity at rotor outlet. 

In addition to these above limitations, there are physical 
limitations imposed by the computer programming itself. These 
are listed in the appropriate appendices. 

2.4 Radial Turbine 
2.4.1 General 

The term radial turbine is used to describe a wide range 
of turbine configurations. In the general sense, it describes 
a turbomachine where the fluid enters the rotor with no appre- 
ciable axial component of velocity. Whether the fluid flows 
radially inward or radially outward and whether it leaves the 
rotor in a radial or axial direction depends upon the context 
of the discussion. The choice made here to limit the radial 
design to that of only one configuration was based on practi- 
cality. The 90 degree inward flow radial turbine, in which 
the fluid enters the rotor in a radial but leaves in an axial 
direction, is the most common type of radial turbine (7). It 
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has the advantage of the strongest structural strength, also 
a factor which can be important at high rotational speeds. 

Only single stage turbines will be considered as these are 
the most common also. 

The methodology of designing the radial turbine can be 
summarized in the following four steps: rotor design, nozzle 
design, scroll design and efficiency correlation. The 
station designations are as shown in Figure 6. 

2.4.2 Output Required 

The parameters necessary to describe the turbine of 
Figure 6 can be narrowed down to the following list for 
preliminary design: 

1) Nozzle Inlet and Outlet diameters 
'2) Nozzle blade height 

3) Rotor Inlet diameter 

4) Rotor Outlet diameters at hub and tip 

5) Maximum Scroll Area 

6) Unit length 

7) Specific speed 

8) Efficiency 

9) Velocity triangles at 

A) Nozzle Inlet 

B) Rotor Inlet and Outlet. 

The velocity triangles for rotor outlet will be described 
at hub, mean, and tip locations. 

2.4.3 Input 

The parameters necessary to describe a radial machine can 
be categorized in the same manner as the axial parameters 
(section 2.3.3): those which describe the mass flow rate, those 
which describe the work accomplished per stuge, and those which 
describe the performance of a stage. 
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Radial Turbine Station Designations 




Machine 



Axi s 



Figure 6 

As it has been previously assumed, that the radial turbine 
design is for a single stage machine, the work accomplished per 
stage is the same as the total work to be accomplished by the 
cycle parameters, section 2.2.3. 

The mass flow rate for this machine can also be described 
in a manner similar to that for the axial machine and which 
is presented in section 2.2.2. The variables necessary for 
this representation are the mass flow rate and the number of 
machines operating in parallel along the fluid flow path. 

Performance estimation is more standardized for radial 
turbines than axial turbines. Through the use of a non-dimen- 
sional specific speed parameter it is possible to correlate 
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efficiency with rotational speed, volumetric flow rate, and 
work accomplished, for turbines designed with the same basic 
assumptions and geometric relationships, (8) There are several 
specific speed parameters in use which will be discussed in a 
later section; however, for the immediate purpose, rotative 
speed is a necessary input. One last input remains and that 
relates to an assumption to be made in the succeeding * section. 
Without explaining this assumption in detail, it is sufficient 
to say that the number of radial blades of the rotor should 
be included as an input. 

The input for the radial design can be summarized as 
follows 2 

1) Mass flow rate of the cycle 

2) Number of units in parallel 

3) Shaft revolutions per minute 

4) Number of radial rotor blades 

5) Cycle conditions of section 2.2. 

2.4.4 Rotor Calculations 

The first step in determining the rotor parameters is to 
determine the conditions at rotor inlet (Station 3). Borrow- 
ing an analysis from the field of centrifugal compressor 
design, it can be shown that the optimum performance is ob- 
tained when the following relationships hold at rotor inlet. (8) 

? = 1 - 2.0 

V 63 = 1 
U„ 



where Z is the number of radial rotor blades. 
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This factor (^) is often called the slip factor. If the 
velocity at rotor exit is assumed to be entirely in the axial 
direction for reasons similar to those given for the axial 
turbine in section 2.3.4, then the slip factor may be used in 
the Euler equation to solve for the radius at rotor inlet. 

Ah = ? U _ 2 

o « 3 

g J 
e o 

where U 3 ='rfRPM r^/lO. 



Substituting the constants into this gives: 



r = 1510.8233 (ah ) 2 . 
J RPM ° 1 



( ? f 



The condition necessary to avoid negative velocities at rotor 
inlet ( 2 ) is: 







2YT. 

Z 



Combini 

possibl 



ng 

e 



this expression with that of 
to determine the optimum inlet 
Tan cKj = Vq 3 = ^ 3/^3 

v v“7U7 

m3 m3 3 

oC = ArcTan ( Z-2 ) . 

3 (2vr) 



the slip factor, 
flow angle, . 



it is 



It is now possible to determine the complete velocity triangles 
at rotor inlet, Figure 7. 

Using the above relationships and the definition for , 



V 



m3 



2jf U 
Z 



U = VfKPM 
J 30 J 

V „ = 994.0819 Uh P 
m3 — — 7} o'jl 

L ( \ ) 2 

V 3 = V m3 /cOS ‘V 

V 63 = V 3 Sin <<K 3 ) 

/ 3> 3 = ArcTan(V 03 -U 3 )/V m3 ) 

V 3 = V m 3 /cos(/? 3 ). 

Radial Turbine Velocity Triangle 



Rotor Inlet 




Rotor Outlet 




Figure 7 
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Using the compressible flow relationships, the flow area 
at inlet is obtained. The flow area and annulus area are 
again related by the cosine of the inlet flow angle as shown 
in Figure 2. The annulus area in this case, however, consists 
of the circumference times the blade height. From this the 
blade height can be calculated. 

At rotor outlet, the assumption is made, as stated pre- 
viously, that the tangential velocity is zero. Watanabe (9) 
states that there is a definite ratio of the mean outlet 
radius to inlet radius which will give the peak efficiency. 
This ratio is approximately 0.6. It has also been somewhat 
traditional to design for a relative outlet angle at the 
shroud of minus 60°. (2) This represents the compromise be- 

tween a low RPM machine and a low loss machine. Applying both 
of these criteria in a modified way, an outlet velocity tri- 
angle at mean, with a relative flow angle of minus 60° was 
chosen. This later proved to correlate well with previous 
work as shown in section 3.6 and avoided an iterative proce- 
dure for calculating rotor outlet conditions. 

Using these assumptions, the outlet velocity at mean, in 
terms of input variables, is 

i. 

V ^ = 9 * 1 347 (* h / 1 ) 2 (Radius Ratio) 

The value of 0.6 specified above is used as a first estimate 
for radius ratio. If the resulting outlet conditions do not 
meet the checks specified in section 2.4.9, then the ratio is 
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increased and outlet conditions recalculated. 

As £ ’ is minus 60°, the other outlet conditions at 

' 4m 

mean are 



U„ =VT RPM r. 

4m 4m 

W. = V. /Cos {£. ). 

4m 4m '4m 

The compressible flow relations are again used to de- 
termine the outlet flow area. In this case, since the outlet 
velocity is in only an axial direction, the flow area is 
equal to the annulus area. The outlet height may be calcula- 
ted by the same method used in section 2.3.5 to calculate 
axial blade height: 

outlet height = Area ann/2y'r^ m . 

Making the assumption that the outlet velocity is con- 
stant with respect to radius, and knowing the outlet hub and 
tip radii, the velocity triangles for these points are calcu- 
lated as previously shown, with the exception that the rela- 
tive outlet angle varies with the radius, or more simply 

ft. = ATAN ( -U./V. ). 

4 4 4m 



2.4.5 Stator Calculations 

The stator calculations are relatively simple. Vatanabe 
(9) states that the clearance between stator outlet and rotor 
inlet should be twice the passage height. if it is assumed 
that the blade or passage height is the same for both rotor 
and stator, then the stator outlet radius is simply the radius 
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at rotor inlet plus two times the blade height. As in most 
cases this gap will not be large, and as it can be assumed 
that the moment of momentum, remains constant across the gap, 




where «sl. 

This means that the velocity triangle at stator outlet is sub- 
stantially the same as at rotor inlet and will not be recalcu- 
lated. 

The radial spacing between blades or pitch, assuming that 
there are as many stator blades as rotor blades, is the circum- 
ference divided by the number of blades. The throat opening 
can also be determined in a manner similiar to that of 
section 2.3.5. Using a pitch to chord ratio of 0.6, the blade 
chord may be calculated. Using the assumption that Cos ^ = 
opening/spacing , the approximate radius at the stator inlet 
can be calculated as shown 

Angle 1 = 

Angle 2 = 

Angle 3 = 

Side 1 = 

Side 2 = 

radius^ = 



in Figure 8. 

hf - Yr 
2 Z 

Jr - oc 
2 

Yr + oC 
Z J 

r~ + C Cos('T[' + ° < ' 

^ z J 

C Sin (If + C<L ) 

Z J 

(Side 1“ + Side 2 “ ) 2 



Radial Turbine Stator Configuration 
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The inlet angle can now be determined using the law of Sines: 

Sin(c^) = r ^ Sin (yr( 1-1,)- ) . 




For the passage from stator inlet to outlet, it can be assumed 
that the fluid is incompressible. This is a valid assumption 
as the Mach number is low at rotor inlet. Since the mass flow 
rate must remain constant, the velocity change is propor tional 
to the flow urea change. This can be expressed as 
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Having determined the inlet angle, c<, the rest of 
triangle is also easily obtainable. 

2.4.6 Scroll Calculations 

The scroll shape assumed is shown in Figure 
assumed that the area is distributed radially in 
relationship, thus : 

Area n — A (1-0 ) 

0 max 2 ^ 



the velocity 



9 . It is 
a linear 



where 0 = radial angle in radians. 

The maximum area needed can be approximated if it is 
assumed that the tangential and meridional components of 
velocity do not vary as a function of the radial angle. The 
maximum area is then equal to the flow area at stator inlet 
which in turn is related by the cosine of the inlet angle. 



Radial Turbine Scroll Configuration 




->Axis of Machine 



Figure 9 
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A =2 r-.Yft cos(c< ) 
max 1 1 

where t is the blade height at rotor inlet 
2.4,7 itadial Efficiency 

As with the axial turbine, the theoretical determination 
of efficiency for a radial turbine is a much debated subject. 
Wood (4) and others have correlated efficiency with the non- 
dimensional specific speed mentioned previously, section 2.4.3 
however, this does not give an exact enough answer for the 
purposes here. A more analytical method is needed. 

Heitt & Johnston (10) relate the performance of several 
specific turbines to several design variables; however, the 
method is rather a narrow one intended to experimentally 
evaluate performance as a function of specific variables, 
rather than provide an analytic efficiency determination. 
Rodgers ( 1 1 ) proposes a method of determining efficiency 
which accounts for losses caused by: 

1) Disc friction 

2) Nozzle shape 

3) Rotor friction 

4) Rotor blade loading 

5) Rotor secondary flow 

6) Rotor incidence 

7) Rotor clearance 

8) Rotor cooling 

9) Exhaust diffuser. 

This method, while verified to within 2 percent, takes into 
account considerably more variables than are necessary for 
the purposes of this mathematical model. 

The method chosen was that used by Rohlik (8) in his 
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study of the effects of turbine geometry on efficiency. This 
method is a compilation of several experimental correlations 
relating the efficiency to: 

1) Stator and boundary layer losses 

2) Clearance losses 

3) Rotor windage losses 

4) Rotor boundary layer losses 

5) Exit velocity losses 

This method provided the needed precision and was based 
on the work of many acknowledged experts in the field; 

Balje (12), Heitt k Johnston (ll), Wood (7), Stewart, Whitney 
& Wang ( 1 3 ) and Shepherd (14) to name a few. The correction 
to account for exit velocity, however, is omitted; thus the 
efficiency calculated is a total to total efficiency and is 
compatible with the efficiency calculated for the axial tur- 
bine. If a total to static efficiency is desired, it can be 
easily calculated from the data given for the outlet velocity 
triangles as: 



7ts - t ]tt 



n + v ^ 

< tt _4m 

2g J^h -,-h 0 ) 

^o ol o2s 

where = efficiency (total to total) 

^ = efficiency (total to static) 

This conversion factor can also be applied to determine total 

to static efficiency for the axial turbine by using in 

lieu of V. . 

4m 
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The mechanics of the method, consist of relating each loss 
to an associated enthalpy change and then summing these in 
a fractional way to determine efficiency. The first loss con- 
sidered is that of rotor windage. Rotor inlet density is re- 
quired and also an absolute viscosity p{) , at rotor inlet. 

The density is computed using the one-dimensional compressible 
flow relationships. 



/O = p 

2 

P = P /RT 
o o' o 

2 i/(r-i)/ 2 4 i/(r-n 

p = P (1- XpL M 2 ) / (1-(<T-1 )" M 4 ). 

w- 2 / 2 

o 



For Mach numbers less than one, the denominator approa- 



ches one and the density is then: 

/° = p (l- d A -i m 2 
/ _ 0 _ — 

RT 



1/ (ar'-i ) 

). 



o 

Using this density, the absolute viscosity, p) can be 
calculated using the kinematic viscosity ("3) given as an 
input . 



A Reynolds number for the rotor inlet is defined as: 

> H = 2U^ r ^ 

e 3 3 * 



Substituting in the constants and the definition of rotative 
velocity: 
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li = 0.2094 395 ltPM r_ /s -■ 
e 3 3 

The windage loss can now be computed as: 

2 



L v = °- 56 /l^l D i 



K ' 2 &. 10 6 

e x 



or 



L = 2,24/y^.r 



w 



3 3 . 



R * 2 m . 10 6 
e x 



The clearance loss is expressed as: 



L = h H 
c o 



where M is defined as 



JA = j ( . 004r^ + .005r 4t ) , 



The rotor boundary layer loss equations involve several 
assumptions (8) which go beyond the scope of this model; 
however, they are consistent with good practice and those 
made previously in this model. It is assumed that the re- 
sulting equations apply to this model. This holds true for the 
stator boundary layer loss equations also. 

The average of inlet and outlet blade spacing is defined 



as : 

S =U(r 3 « |h +r 4t ). 



The solidity is defined as: 



<3"= 1.6r ((r /r )-l). 

4 m 4 m 3 
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A loss coefficient ^ e T ) is calculated by: 

e = ( 0.017 c~) 



(1 + 1. 9S ) 

(0.5-0. 003Z-0 .017 ) (t 3 +t 4 ) 

and the total boundary layer loss is given by: 

2 



L = e 
r r 



W 



4m 



1-e 2g J 
r fe o 

In a like manner for the rotor, a stagger angle is defined 
as : 

where is given by 
o<?~ = Ar cTan ( Sinc^ ) 

( 2t 3 + Cos ' x 3) 



The loss coefficient (e ) is calculated from 

s 

e = 0.0076 (1+Coscc ) 

S Cos- 3-0.025 ( -oTT 1 )’ 

and a total loss for the stator boundary layer calculated 
from 



L = e 
s s 



(v 3 2 ) 

l 7 ^ (2 ^o J) 



The efficiency is found from 



v> .. = ^h -L -L A^h +L +L )• 
7 tt o w c' o s r 

2.4,8 Miscellaneous Parameters 



Specific speed seems to be defined differently by each 
researcher. Some definitions are non-dimensional and others 
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are not. Wood (7) and Kohlik (8) define their specific speed 
similiar to that used in centrifugal pump theory as 

n s = RPM(g 4 )* 

H 3 '' 4 

where is the volumetric flow/sec 

H is the ideal head drop across 
the turbine in feet. 

This, however, is dimensional and for that reason the 

non-dimensional specific speed recommended by Carmichael (2) 

is used instead. 

This is defined as 

N = RPM Q * 

S 60 _5 

( g J Ah ) 

' fo o o 

The rotor outlet volumetric flow can be calculated as 



3/4 



where is obtained using the compressible 
flow relationship of the previous section. 

As a matter of interest the specific speed of Wood et al 

may be obtained from that used in the model by using the 

following formula: 

N (V»ood) = 810.04455 N (Carmichael), 
s s 

The axial unit length of the rotor is a parameter which 



must be calculated in detail for the final design; however, in 
order to simplify the calculations in the mathematical model 
a simpler method is used. Based on the data given in several 
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of the references, on the average, the unit axial length is 
on the order of the inlet radius less the outlet hub radius, 
thus for the purposes here, 

length = r 3 -r 4h . 

The machine diameters given as output are simply calcu- 
lated from the previously derived radii. 

2.4.9 Radial Limitations 

The specific speed parameter has an important influence 
on the shape of the turbine. In order that the model not 
generate dimensions which were unacceptable from a geometric 
point of view (i.e. a radical unconventional design in which 
other factors such as shroud clearance might have an unaccount- 
ed for effect), the input parameters are first checked to de- 
termine the specific speed of the proposed turbine. If this 
specific speed is below 0.05 or above 0.25 the input para- 
meters will not produce a design which is conventionally 
feasible. 

The calculation of stress in the rotor is a complicated 
subject and one which is beyond the scope of the preliminary 
design phase. Wood (7), however, provides a means of deter- 
mining whether a design will be stress limited. Specific 
speed is written in terms of outlet area, RPM, rotor tip speed 
and several fluid velocities. Maximum expected values for a 
stress limited case are then used to solve for a specific 
speed. This specific speed was approximately twice that 
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which would be allowed by the previous criteria and, therefore, 
the assumption is made that the radial turbines allowed by 
this model will not be stress limited. 

As in the axial turbine a continuity check must be made 
at rotor exit to assure that negative diameters are not calcu- 
lated. The rotor outlet radius is checked to insure that it 
is at least 10 percent of the inlet radius. This differs 
from the previous 50 percent used in the axial turbine; however, 
the radial rotor doesn't have to incorporate a through shaft 
and thus can be of a much smaller inner diameter. If the 
radius in question is too small, the radius ratio used in 
calculating V , section 2.4.4, is increased by 2 percent. 

This has the effect of increasing the density of the fluid 
and thus decreasing the outlet annulus area required. 

The calculations described here for the radial turbine 
do not account for any shock effects caused by operation at 
a Mach number greater than 1. liven without this correlation, 
however, Rodgers (11) states that Mach number effects are 
not important until >1 = 1.2. Carmichael (2) recommends a 
limit of 1.1. The lower limit was used to be consistent 
with other assumptions made previously. The two velocities 
likely to be the largest, the outlet tip relative velocity, 
and the rotor inlet velocity are checked to insure that they 
are below a Mach number of 1.1. If not, the input parameters 
result in an unfeasible turbine. 
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EVALUATION OF THE MODEL 

3.1 General 

In order to facilitate- the preliminary design of turbines 
using the mathematical model developed in the previous chapter, 
it was necessary first to program the model, and then to run 
the program over a wide range of inputs in order to verify the 
results with existing correlations; or in the case of the sim- 
pler geometric relationships, with hand calculations. The 
areas which were examined for validity and the methods used 
to establish same can be summarized as: 



Area 

Simple mathematical 
relationships 

Aspect ratio selection 



Soderberg correlation 
function 

Axial turbine efficiency 

f 

Radial turbine efficiency 



Method 

hand calculations 



hand calculations of 
existing turbine 

comparison with 
original function 

comparison with other 
published correlations 

comparison with other 
published correlations 



Each of these areas will be described in more detail in the 
following sections. 

The cycle chosen for the validation was the Zener sea 
cycle as proposed by Clarence Zener (15). This cycle is 
similiar to others (16) which seek to take advantage of the 
temperature differences existing in the world’s oceans, and 
is representative of the type of cycle for which the mathe- 
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matical model was designed to be used to evaluate. The de- 
tails of this cycle and of the individual parameters which 
were examined are contained, in Appendix F. It is sufficient 
to say here that over 50 different combinations of input para 
meters were examined for each of two different working fluids 
This range of input variables provided a representative 
picture of the design surface, and allowed the results to be 
broad enough to correlate with existing data over a wide 
range , 

3.2 Simple Mathematical Relationships 

This facet of the model is the predominent one, as 
explained in section 1,3, and also the easiest to verify. 

The major portion of the model is based on traditional theory 
which has been well proven in the laboratory. It is not the 
purpose here to attempt to revalidate the conclusions of re- 
searchers far more qualified than the author. The emphasis 
has been placed instead on choosing accepted relationships 
and correlating them into a usable model. 

Based on this, the validation of this segment consisted 
of hand calculating the various mathematical relationships 
and verifying the results of the computer program. To this 
end, the model has been debugged of all forseeable mechanical 
type mathematical errors, 

3.3 Aspect Ratio Correlation 

The decision to develop an internal choice of aspect 
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ratio based on blade bending stress for the axial turbine, 
rather than have the designer predetermine it, was based on 
primary and secondary input, factors as explained in section 
2.3.3. It was necessary then to verify that the internal 
choice was a reasonable one, consistent with existing designs. 

This was carried out by applying the correlation developed 
in section 2.3.4 to an existing design for an intermediate 
stage turbine, developed for marine propulsion use and des- 
cribed by T. B. Hutchinson in reference 17. This design was 
for an 8 stage steam turbine operating at 5036 RPM and pro- 
ducing 8940 horsepower at 90 percent efficiency with an inlet 
temperature and pressure of 1000°F and 225 psig respectively. 
Using the data supplied by Hutchinson in reference 17, and 
making assumptions, where necessary, of input variables which 
would optimize efficiency, the correlation of section 2.3.4 
was applied to the design and an aspect ratio was computed. 

This aspect ratio was approximately 1.4. The aspect ratio 
chosen in the actual design was approximately 1.5. 

While this simple point verification does not guarantee 
validity over the complete range of input variables, it does 
show that the relationship is reasonable. When this relation- 
ship is combined with an upper limit for aspect ratio, to take 
into account the fact that other factors than stress may limit 
the design, the resulting range of variation of aspect ratio 
(0.25 to 10. 0) is consistent with existing designs (5) and is 
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assumed to provide a reasonable result. 

3.4 Soderberg Correlation 

In determining efficiency for the axial turbine, section 
2.3.5, it was necessary to formulate a function which approx- 
imated the curve of emperical data given by Soderberg (4) for 
the correlation of gas deflection angle and the primary loss 
factor. This was accomplished in a computer subprogram 
through the use of LaUrange polynomials to curve fit the 
given data to a third order approximation of the function 
( 1 8 ) . Without explaining the details of the theory of curve 
fitting, it is sufficient to state that within a given sub- 
interval the correlation was represented by the function 
W = X x Y 1 + X 2 Y 2 + X 3 Y 3 

where X^ , X^ and X 3 are third order functions dependent upon 
the value of the abscissa given ? and Y^ , Y^ , and Y 3 are values 
of the ordinate which bracket the given abscissa value. 

The simulation of this function is shown' in Figure 10. 

3.5 Axial Turbine Efficiency 

Although the correlation adopted for determination of 
axial total to total efficiency has been previously validated, 
it was necessary to compare the results obtained with the 
model to other correlations in order to assure that the 
correlations provided reasonable results for the assumptions 
made in the model. This was accomplished by comparing the 
computed efficiency of several example cases for the Zener 
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Subprogram SODCOR Simulation of Soderberg ^ Parameter 




Gas Deflection in Degrees 



I 



* 



Figure 10 
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cycle with the correlations presented by Craig and Cox (3). 
This comparison is shown in Figure 11. 

This comparison does not account for the differences 
between the two correlations in their accounting for secon- 
dary loss parameters; however, the results are close enough 
to show that the application of Soderberg T s method is valid 
and consistent with accepted practices. 

3.6 Radial Turbine Efficiency 

It was also necessary to compare the computed radial 
turbine efficiency to other correlations in order to verify 
its application in the model. This can be more easily accom- 
plished in the radial turbine than the axial through the use 
of the non-dimensional specific speed parameter mentioned 
in section 3.4.8. Several examples, covering a wide range 
of specific speeds, were compared with the correlations of 
Wood (2) and Iiohlik (8). The results of this correlation 
are shown in Figure 12. The close agreement over the wide 
range of specific speeds shows that the determination of 
radial efficiency is also valid and consistent with accepted 
practice . 
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Axial Turbine Efficiency Comparison 



Craig & Cox 



Correlation 



Program Output 




Work Coefficient 

) 



j 



Figure 11 
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Radial Turbine Efficiency Comparison 
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CONCLUSIONS 

The formulation of conclusions about a project of this 
nature is a somewhat limited undertaking. The object of the 
thesis is to compile existing data, relationships, and methods 
in the field of turbine design and to formulate and compute- 
rize a mathematical model, based on this compilation. This 
has been accomplished as shown in chapter 2. 

The validity of the mathematical model has been estab- 
lished in chapter 3 to the extent that it is valid, consistent 
with the assumptions stated in its formulation, and within 
the requirements posed by the preliminary design stage. 

The actual use of this programmed model, as explained 
in Appendix F, has shown that through the utilization of 
this program, or programs of this type, the designer may 
explore a wide range of input variables in a comparably 
short amount of time. This enables him to more fully examine 
the design spectrum and to develop a qualitatively better 
estimation of thermal cycle feasibility. 
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RECOMMENDATIONS 

In developing the mathematical model described in chapter 
2, it became apparent that there were four areas where the 
model could be expanded which would increase its usefulness 
to the designer. The first of these is in the area of the 
flow relationships utilized. The one dimensional flow rela- 
tionships do not completely describe the flow. Expansion of 
these relationships to include shock effects would broaden 
the applicability of the model. 

The second and third areas which could be expanded are 
integrally related to one another. They entail enlarging 
the areas of the machine which are described by the model, 
perhaps including material limits and selection, disc design, 
bearing and shafting design. Concurrent with this would be 
the increasing number of input variables necessary to des- 
cribe the machine. This would be limited by the amount of 
detail required to satisfy a particular designer's definition 
of preliminary design. 

The last area which was not included in this model is 
that of cost estimation. Information in this area is gene- 
rally proprietory and not easily obtainable. Most of the 
information available is limited in scope and not applicable 
to the general case as described herein. This area in itself 
could be the subject of a dissertation, but certainly more 
work is needed before it is acceptable for inclusion in a 
model similar to that developed in chapter 2. 
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APPENDIX A 



A.l Variable List 

The following is a list of variables as they appear in 

the main program and the subprograms. They are listed in the 

order in which they appear within the body of the programs. 

Subscripted variables refer to machine inlet and outlet, 

subscripts 1 and 2 respectively. The units of these variables 

are those listed in the Notation section, pages 7 to 9, with 

the exception of pressure, kinematic viscosity and all angles 

2 2 

which are internally converted to lbf/ft , ft“/sec and 
radians respectively. 

MAIN PROGRAM 

Cycle and Machine Inputs 

inlet stagnation enthalpy 
stagnation temperature 
stagnation pressure 
kinematic viscosity 
outlet stagnation enthalpy 
outlet isentropic stagnation enthalpy 
outlet isentropic stagnation temperature 
outlet isentropic stagnation pressure 
specific heat at constant pressure 
ratio of specific heats 
specific heat at constant volume 



HOI 

TO 

PO 

NU 

1102 

H02S 

T02S 

P02S 

CP 

GAMA 

CV 
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H 


universal gas constant 


DHO 


stagnation enthalpy drop across machine 


RJ 


universal gas constant times 778 


MM 


number of combinations of input parameters entered 


AXN 


number of turbine stages-axial 


REAM 


reaction at mean 


RPMX 


shaft revolutions per rninut e-axi al 


ALFAMN 


rotor inlet true flow angle at moan 


U'NITX 


number of axial turbines in parallel on flow path 


RPMR 


shaft revolutions per minute-radial 


UN1TR 


number of radial turbines in parallel on flow path 


ZR 


number of radial blades 


WW 


total mass flow rate of working fluid 



Axial Turbine 



itADM 


radius at mean 


VTMN 


tangential velocity component at rotor inlet, 
mean radius 


VwN 


true velocity at rotor inlet, mean radius 


VX 


axial velocity component 


Wj'iR 


relative velocity at rotor outlet, mean radius 


TALFMN 


tangent of ALFAMN 


BETAMN 


relative flow angle at rotor inlet, mean radius 


BETAMit 

ZETAN1 


relative flow angle at rotor outlet, mean radius 

primary Soderberg efficiency parameter at rotor 
inlet 



ZETAN1 



ZETARi 


primary Soderberg efficiency parameter at rotor 
outlet 


CON 


product or radius times tangential velocity 
component 


WMN 


relative velocity at rotor inlet, mean radius 


UMN 


tangential rotor velocity at mean radius 


ELEN 


function f (section 2*2.3) 


AELO 


flow area 


AANN 


annulus area 


T 


blade height-axial 


RADT 


tip radius 


RADII 


hub radius 


VTTN 


tangential velocity component at rotor inlet, 
tip radius 


VTHN 


tangential velocity component at rotor inlet, 
hub radius 


ALEATN 


rotor inlet true flow angle at tip radius 


ALEAHN 


rotor inlet true flow angle at hub radius 


VTN 


true velocity at rotor inlet, tip radius 


VHN 


true velocity at rotor inlet, hub radius 


UTN 


tangential rotor velocity at tip radius 


UHN 


tangential rotor velocity at hub radius 


WTTN 


relative tangential velocity component at rotor 
inlet, tip radius 


W TUN ' 


relative tangential velocity component at rotor 
inlet, hub radius 


RETATN 


relative flow angle at rotor inlet, tip radius 


BETAHN 


relative flow angle at rotor inlet, hub radius 
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WTN 


relative velocity at rotor inlet, tip radius 


WHN 


relative velocity at rotor inlet, hub radius 


BETATR 


relative flow angle at rotor outlet, tip 


radius 


BETAHR 


relative flow angle at rotor outlet, hub 


radius 


WTR 


relative velocity at rotor outlet, tip radius 


Willi 


relative velocity at rotor outlet, hub radius 


CHI) 


blade chord 




ASliAT 


aspect ratio 




SPN 


stator spacing 




TBTAMN 


tangent of BETAMN 




TBTAMR 


tangent of BETAiiR 




SPR 


rotor spacing 




TR 


throat opening for the rotor 




TN 


throat opening for the stator 




DR 


hydraulic diameter for rotor 




DN 


hydraulic diameter for stator 




RER 


Reynolds number for rotor 




REN 


Reynolds number for stator 




ZETAR2 


Soderberg efficiency parameter corrected 
aspect ratio for rotor 


for 


ZETAN2 


Soderberg efficiency parameter corrected 
aspect ratio for stator 


for 


ZETAR3 


Soderberg efficiency parameter corrected 
and Re, for rotor 


for AR 


ZETAN3 


Soderberg efficiency parameter corrected 
and Re, for stator 


for AR 


EPP 


stage efficiency 
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ZX 

NZX 

EFFX 

XMN 

XMR 

UEAH2 
ESTMX 
DI AMI 
DJLAM2 
UNLExM 
(i 

itadial 

Z 

HAD 3 

ALFA3 

VM3 

V3 

XM3 

FUN 3 

AFL03 

Tit 

KAU2 

itiiAT 



number blades-decimal number 
number blades-whole number 
machine efficiency 

true Mach number at rotor inlet, hub radius, last 
stage 

relative Mach number at rotor outlet, tip radius, 
last stage 

reaction at last stage, hub radius 
centrifugal stress at blade root, last stage 
inlet diameter 
outlet diameter 
machine length 

radians to degrees conversion 

Turbi ne 

slip factor 

radius at rotor inlet 

true inlet flow angle, rotor inlet 

meridional velocity component at rotor inlet 

velocity at rotor inlet 

Mach number at rotor inlet 

function f at rotor inlet (section 2.2.3) 

flow area at rotor inlet 

stator blade height 

radius at stator exit 

radius ratio (section 2.4.4) 
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V 4i v i 


velocity at rotor outlet, mean radius 


V4H 


velocity at rotor outlet, hub radius 


KAD4M 


mean radius at rotor outlet 


XM4M 


Mach number at rotor outlet 


FUN 4 


function f at rotor outlet (section 2.2.3) 


AANN4 


annulus area at rotor outlet 


T4 


blade height at rotor outlet 


RAD4H 


hub radius at rotor outlet 


RAD4T 


tip radius at rotor outlet 


BB 


exponential constant 


BEN4 


fluid static density at rotor outlet 


Q 


volumetric flow rate at rotor outlet 


SPSP 


non-dimensional specific speed 


SPRR 


stator outlet spacing 


OP 


stator throat opening 


UR 


stator chord 


RADI 


stator inlet radius 


SALFA1 


sine of ALFA1 


UALFA1 


cosine of ALFAl 


ASCR 


maximum scroll area 


VI 


velocity at stator inlet 


ViMi 


meridional velocity component at statoi' inlet 


V'i’l 


tangential 


ALFA1 


true flow angle at stator inlet 


VT3 


tangential velocity component at rotor inlet 
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U3 


tangential rotor velocity at inlet 


BETA3 


relative flow angle at rotor inlet 


V3 


relative velocity at rotor inlet 


U4M 


tangential rotor velocity at outlet, mean radius 


U4H 


tangential rotor velocity at outlet, hub radius 


U4T 


tangential rotor velocity at outlet, tip radius 


BETA4M 


relative flow angle at rotor outlet, mean radius 


BETA4H 


relative flow angle at rotor outlet, hub radius 


BETA4T 


relative flow angle at rotor outlet, tip radius 


W4M 


relative velocity at rotor outlet, mean radius 


W4H 


relative velocity at rotor outlet, hub radius 


W4T 


relative velocity at rotor outlet, tip radius 


XM4 


relative Mach number at rotor outlet, tip radius 


UE3 


Reynolds number at rotor inlet 


DEN3 


fluid static density at rotor inlet 


XLM 


windage loss 


XLC 


clearance loss 


SPltM 


mean rotor spacing 


SiGMAR 


sol idi ty 


Eli 


rotor boundary layer loss coefficient 


XLR 


rotor boundary layer loss 


ALFAO 


stator angle (section 2.4.7) 


ALFAST 


stator stagger angle 


ES 


stator boundary layer loss coefficient 


XLS 


stator boundary layer loss 
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EFFIt 


machine ef f iciency 


UNLENli 


machine length 


DIA3 


diameter at rotor inlet 


DIA4T 


tip diameter at rotor outlet 


DIA4M 


hub diameter at rotor outlet 


D1A1 


diameter at stator inlet 


DIA4M 


mean diameter at rotor outlet 


DIA2 


diameter at stator outlet 


SUBPROGRAM 


SODCOR 


XI 


gas deflection angle, radians 


X 


gas deflection angle, degrees 


XX 


abscissa for Soderberg correlation 


II 


ordinate 


XL1 


first LaGrange Polynomial 


XL 2 


second LaGrange Polynomial 


XL 3 


third LaGrange Polynomial 


I 


Soderberg gas deflection loss parameter 


SUBPROGRAM 


FLOREL 


GAMA 


ratio of specific heats 


R 


universal gas constant 


CP 


specific heat at constant pressure 


TO 


stagnation temperature 


V 


velocity 



V 



xm 

FLFN 

T 

EXPll 



(JON 



Mach number 
function f (section 
static temperature 
internal constant 
internal constant 



FUN 



internal constant 
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APPENDIX B 



B.l i v lain Program 



A complete listing of .the main program is contained on 
the following pages. Appendix E contains detailed instruc- 
tions for input-output and running the program. The mathe- 
matical model for the program is developed in chapter 2. 

The logic sequence used in programming the model is outlined 



1 ) 

2 ) 

3 ) 

4 ) 

5 ) 

6 ) 

7 ) 

8 ) 
9 ) 

10 ) 

11 ) 

12 ) 

13 ) 

14 ) 

15 ) 

16 ) 

17 ) 

18 ) 



19 ) 

20 ) 
21 ) 



Definition of real variables and common variables 

Dimensioning of subscripted variables 

Input format specification 

Output format specification 

Keadin input data 

Readout input data for checking 

Convert input data to standard units 

Calculate fluid properties 

Readout fluid properties for checking 

Commence do-loop for calculations on each set of 

turbine input parameters 

Calculate axial turbine mean parameters 
Calculate axial first and last stage parameters on 
internal do-loop 

If design is rejected on specified checks, axial cal 
culations are dropped and radial are commenced. 
Calculate axial machine parameters 
Convert axial variables to output units 
Readout output data for axial machine 
Calculate radial parameters 

If design is rejected on specified checks, radial 
calculations are dropped and next set of axial 
calculations coinmenc ed 

Convert radial variables to output units 
Readout output data for radial machine 
Commence calculations for next set of axial input 
parameters 
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APPENDIX C 



C Subprogram SODCOR 

This program computes 'an efficiency parameter from 
given data points using LaGrange Polynomials in a curve- 
fitting logic. A complete listing of the subprogram is 
contained on the following page. The mathematical model 
is explained in section 3.4. The logic sequence used in 
programming the model is outlined below: 

1) Transfer of input variables 

2) Definition of curve data points 

3) Conversion of input data to standard units 

4) Classification of input data as to sub-division 
of curve for curve-fitting purposes 

5) Definition of sub-division curve-fitting data 

6) Calculation of LaGrange Polynomials 

7) Calculation of output variable 

8) Transfer of output variable 



S|'|BR0UTINE SoDCOR (Xi,Y) 
Dimension yX( 1 * >,yy( is ) 
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APPENDIX D 

D Subprogram PLOUEL 

This subprogram computes a flow function and Mach number 
using the one-dimensional flow functions. The mathematical 
model is explained in section 2.2.3. A complete listing of 
the program is contained on the following page. The logic 
sequence used in programming the model is outlined below. 

1) Transfer of input variables 

2) Transfer of common variables 

3) Calculation of Mach number 

4) Calculation of flow function 

5) Transfer of output variables 



96 



z 


















U 










— * 








3 










51 








Li. 






— 




X 














* 


**» 


& 








21 


•*-« 




(V 


a: 


J~ 


— 






X 






1 


x 


X 


2 








CL 




<- 


IX 


* 


O 






> 


U OJ 




n 


-4- 


— 


U 






•v 


* ix 




c 


LD 


• 


-fr 






0 


-4 CO 




0 


CO 


CM 


Z 






h- 


■4- • 




* 


*-» 


X 


X 






— - 


ix *-< 




• 


•3“ 


- — ■ 










• co 




OJ 


G. 


• 


X 






UJ CL 


OJ cs 




— - 


e 










a; 0 


xC LT, 




X 


Cl. 


t 


kJ 






0 <* 


IS OJ 





— ■ 


* 


< 


£1 






-J Ct 


ts * 


*— 


• 


<i 




X 






Ll 


LO L- 


3 


«— i 




< 


UJ 






<x 


— * 


O 


+ 


< 




* 






Lu r 


\ Cr 


— 


<a 


rr* 


— - 


* 






2 < 


> * 


I— 








z 






►h e; 


* c 


CL 


< 


L— 




3 






*— 


> 2. 


O 


O 


LX 


+ 


Ll 






3 z 


— < 


OJ 




O 


c* 




-> 




0 0 


• O 


X 


<1 


(/) 




II 


a: 




cr >: 


O H 


> 


GC 


II 


1) 


2* 






CG 5 : 




a 


ll 


21 


21 


u. 


i — 


0 


“ c 


'11-3 


2 


X 


c- 


■ _ ; 


: ^ 1 


► Ll 


z 


CO u 


1“ O 


X 


UJ 


u 


Ll ' 


Ll 


cr 


UJ 



97 



APPENDIX E 

E.l Specific Instructions for Use 

This section is intended to provide detailed instruc- 
tions for the use of the computer programs described in the 
previous appendices including input-output interpretation. 
The specific program control cards, described are required 
on the Interdata 70 computer. Check the user’s manual be- 
fore running on other machines as these control cards may 
differ . 

K.2 Loading Sequence 

The sequence for loading the main and subprograms with 
their required control cards is as follows: 

1) Job card 

2) Language card 

3) Listing designation card 

4) Main program 

5) Subprogram SODCOR 

6) Subprogram PLOREL 

7) Execution card 

8) Common length designation card 

9) Data cards 

10) Termination card 

E.3 Job Control Cards 

* 

' The JOB card, language card, listing designation card 
and JOB ’ terminati on cards are all of the standard format, 
common to any job run on the Interdata 70 system. The ex- 
ecution card must carry a numeral one in column seventeen in 
addition to the standard format. The common length designa- 
tion card specifies the space in blank common to be used in 
the job. The format for tiiis card is *BC 000C commencing 
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in column one, 

E. 4 Data Cards 
E . 4 . 1 General 

The Data cards necessary for use of the program can be 
grouped into three categories. 

1) Cycle Data 

2) Number of Design Sets 

3) Turbine Design input sets 

They must be loaded in the above order. Card groups 1 and 2 
are only required once per job, while card group 3 is re- 
quired for each individual design set examined. 

E.4.2 Cycle Data 

Three cards are required to specify the data for a 
given thermal cycle. The context of each is given below. 

No information can be omitted, and they must be loaded in 

# 

the order described. 



SPACES 



INFORMATION 



UNITS 



FORMAT 



CARD ONE 



1-10 

20-30 

40-50 

60-70 



inlet total enthalpy 
inlet total temperature 
inlet total pressure 
inlet kinematic viscosity 



BTU/lbm 

°R 

psia 



F10.5 

F10.5 

F10.5 

F10.5 




CARD TWO 



1-10 

20-30 

40-50 

60-70 



outlet total enthalpy BTU/lbm 

outlet total temperature °R 

outlet total pressure psia 

2 

outlet kinematic viscosity ft /hr 



F10 . 5 
F10.5 
F10.5 
FI 0.5 
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SPACES 



INFORMATION 



UNITS 



FORMAT 



CARD THREE 



1-10 



outlet isentropic total BTU/lbm 
enthalpy 



F10.5 



20-30 



outlet isentropic total °R 

temperature 



F10.5 



40-50 



outlet isentropic total 
pressure 



psia 



F10.5 



E.4.3 Number of Design Sets 

One card is necessary to specify the number of sets of 
turbine parameters entered. A set consists of one combina- 
tion of axial turbine input parameters, and one combination 
of radial turbine input parameters, with both combinations 
sharing the same total mass flow rate (section 2.2.2). The 
program will accommodate up to 99 separate sets of input 
data. The number of sets inusf be specified in 12 format on 
spaces one and two of the card with the other spaces left 
blank. This card may not be omitted. 

E.4.4 Turbine Design Input Sets 

A turbine design input set consists of a combination of 
axial input parameters, section 2.3.3, and a combination of 
radial input parameters, section 2.4.3, which share the same 
total mass flow rate, section 2.2.2. Three cards are 
necessary for each set and no information or card may be 
omitted. They must be loaded in the order in which they are 
described. in the event that a run is not to be designated 
by name, a blank card must be included in place of card one. 
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SPACES 


INFORMATION 


UNITS 


FORMAT 


CARD ONE 








6-70 


Any desired run designa- 
tion information 




H 


CARD TW'O-Axial Turbine 






11-15 


Number of stages 




F5.0 


21-25 


Inlet reaction at raidblade 




F5.2 


31-40 


Shaft RPM 


Rev/ rain 


F10.0 


51-55 


Rotor inlet flow angle 
at midblade 


Degrees 


F5.2 


61-65 


Number of machines op- 
erating in parallel on 
total flow path (may be one 


) 


F5.0 


CARD THREE- 


-Radial inflow Turbine 






11-20 


Shaft RPM 


Rev/ min 


F10.0 


31-35 


Number of machines op- 
erating in parallel on 
total flow path (may be one 


) 


F5.0 


41-45 


Number of rotor blades 




F5.0 


51-60 


Total mass flow rate for 
all machine s-i . e . sum of 
all parallel flow paths. 


lbra/ sec 


F10.2 



See section 2.2.2, It 
must be the same for Axial 
and Radial Design Set. 

E.5 Output 
E , 5 . 1 General 

The cycle data will be printed out once for each job, 
in the same format it was entered in with, as a check for 
accuracy. The universal gas constant, in units of ft Ibf/ 
lbm°R, the specific heat at constant pressure, and the ratio 
of specific heats will also be printed out as a check. For 
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each turbine design data set, the run designation information, 
the assumed input variables, and the preliminary design data 
will be printed out. If the assumed input variables result 
in a design which doesn’t pass the axial or radial feasibi- 
lity checks, sections 2.3.8 and 2.4.9, the design data will 
be omitted and an unfeasibility message printed, section 
E.5.3. 

E.5.2 Preliminary Design Data 

The format of the preliminary design data is self- 
explanatory. The output is in units of feet for all dimen- 
sions, degrees for all angles, and feet per second for all 
velocities. Stress is given in psi and tpi . The abbrevia- 
tions used in describing the velocity triangles are as 
follows : 

Axial Turbine 

UHN, UMN, UTN, - rotor tangential velocity at hub, mean and 

tip radii 

VTHN , VTMN, VTTN, - fluid tangential velocity component at 









rotor 


inlet at 


hub 


> 


mean 


and tip 


radi i 


vx 






- fluid 


axial ve 


loci 


ty 


component 




VHN, 


VMN, 


VTN, 


- fluid 


velocity 


at 


ro 


tor 


inlet at 


hub, mean 








and ti 


p radii 












WHN, 


WiMN , 


WTN , 


- fluid 


relative 


vel 


oc 


i ty 


at rotor 


inlet at 








hub, mean and 


ti p 


ra 


di i 






Willi , 


>lMli , 


WTli , 


- fluid 


r e 1 a t i v e 


vel 


oc 


i ty 


at rotor 


o utl e t , hub , 



mean and tip radii 

ALEAI1N , AEPAMN , ALFATN - rotor inlet flow angle at hub, mean and 

tip radii 
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BETAHN,BETAMN,BETATN, - relative rotor inlet flow angle at 

hub, mean and tip radii 

BETAHR ,BETAMU , BETATR , - relative rotor outlet flow angle at 

hub, mean and tip radii 



Radial Turbine 



VI, V3 


— - 


fluid velocity at station 1 or 3 


VT1, VT3, 


- 


fluid tangential velocity component at 
station 1 or 3 


VM1, VM3, 


- 


fluid meridional velocity component at 
station 1 or 3 


U3 


- 


rotor tangential velocity at inlet 


V3 


- 


fluid relative velocity at rotor inlet 


ALFA1 , ALFA3 , 




inlet flow angle at station 1 or 3 


BETA3 




relative rotor inlet flow angle 


U4H,U4M,U4T, 


- 


rotor outlet tangential velocity at hub, 
mean and tip radii 


V4H, V4M,V4T, 


- 


rotor outlet fluid velocity at hub, mean 
and tip radii 


V4H,W4x4,W4T, 




rotor outlet fluid relative velocity at 
hub, mean and tip radii 



BETA4H ,BETA4M ,BETA4T , - rotor outlet fluid relative flow 

angle at hub, mean and tip radii 

E.5.3 Unfeasibility Conditions 

If a given set of turbine design inputs does not meet 

the feasibility checks described in sections 2.3.8 and 2.4.9, 

a design unfeasibile message will be printed. As a reference 

the checks associated with each message are summarized below. 

MESSAGE INTERPRETATION 

Axial Turbine 



* 



Either ^ or syS, ) is greater 
than 140° 



Gas Deflection too big 
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MESSAGE 


INTERPRETATION 


Hub Radius too small 


r h <0 - 5 r t 


Hub iieal Mach Number 
too big 


Last stage hub true Mach 
number greater than 0.99 


Tip Real Mach.Nutnuer 
too big 


Last stage tip relative Mach 
number greater than 0.99 


Hub reaction too small 
Radial Turbine 


Last stage hub reaction less 
than 0.01 


Specific Speed to Radical 


Machine specific speed less 
than 0.05 or greater than 

0.25. 


Inlet Mach number too 
high 


Mach number at rotor inlet 
greater than 1.1. 


Outlet Relative Mach 
number too high 

E.5.4 Compilation Errors 


Relative Mach number at rotor 
tip greater than 1.1 



It is possible, given certain sets of axial turbine in- 
put data, to cause a computer run time error (negative square 
root and overflow). This condition sometimes occurs during 
the calculation of aspect ratio. It lias been forseen in pro 
gramming the model and will be internally corrected if it 
occurs. The data output is still valid if this error occurs 
even though an error message is printed out. Other run time 
errors or errors in computing the radial turbine data should 
be cause for re-examination of the input data provided. 
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APPENDIX P 

EVALUATION OF THE ZENER CYCLE 

F.l Coneral 

The Zener Cycle, as described by Clarence Zener (15), 
is based upon the thermocline existing in the oceans of the 
world. According to Zener, it should be possible to take 
advantage of this temperature difference, (20°C is used as 
an example) to power a thermal cycle which could in turn 
produce useful energy. This supposition is based upon the 
assumption that technology could produce the necessary 
equipment. Assuming that this is possible, Zener estimates 
that n the tropical oceans in the year 2000 could supply the 
whole world with energy at a per capita rate of consumption 
equal to the U.S. per capita rate in 1970, and suffer only a 
one-degree C drop in temperature”. This is estimated to be 
about 60 billion kilowatts. Zener proposes a thermal cycle 
using ammonia as the working fluid and operating between 5°C 
and 25°C. His estimated plant size is about 8000 cubic feet. 

The cycle consists of taking high pressure liquid at 
25°C, expanding it through a turbine, and condensing the 
vapor at 5°U. The fluid would then be pressurized and 
heated to 25°C to begin the cycle again. Sea water would 
be used to both heat and cool the fluid through the use of 
heat exchangers. Ammonia is suggested as a working medium. 
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F.2 Working Fluid 

Without exploring the heat exchanger problems involved 
in detail, the feasibility of the turbines needed for such 
a cycle was evaluated. The calculations were based upon a 
plant output of approximately 100,000 kilowatts with an 
assumed turbine efficiency of 90 percent. Based upon this 
the flow rate for a specific fluid was calculated. Ammonia 
and Freon 21 were chosen as the fluids; ammonia because it 
was proposed by Zener, and Freon 21 because it is typical of 
"the recently developed refrigerating fluids" mentioned by 
Zener. The cycle operating points are summarized as follows, 



along with the necessary mass flow rates. 



Point 


Enthalpy 


Temp 


Pressu 




(BTU/lbm) 


(°R) 


(psia 






AMMONIA 




1 


628.74 


528 


124.4 


2 


612.575 


510 


89.19 


2S 


610.779 


510 


89.19 






FREON 21 




1 


127.56 


528 


22.2 


2 


124.323 


510 


15.33 


2S 


123.963 


510 


15.33 






MASS FLOW RATE 






Ammonia - 6050 lbm/sec 





Freon 21 - 30220 lbin/ sec 
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P.3 Range of Input 

The inputs chosen for the radial turbine were based 
upon a back calculation from specific speed to RPM and number 
of units in parallel. The RPM was quantified as being 3600/N 
where N is an integer number. This was done because it is 
assumed that the machine would be used to produce 60 hertz 
power without benefit of a reduction gear. The number of 
units was then specified to provide a variety of machines 
covering the range of specific speed recognized as feasible, 
section 2.4.9. Twelve rotor blades were specified for the 
ammonia trials and 15 for Freon 21. This brackets the usual 
design range. 

The axial turbine inputs were based on combinations of 
parameters used to achieve high efficiency. A large slow 
machine was favored with high reactions and few stages. 

This was not immediately apparent at the outset, but after 
several runs with no combinations yielding a feasible turbine, 
it was indicated. 

A summary of design inputs for one run is shown below. 

The feasibility column indicates whether the design meets 
the feasibility criteria of the model, sections 2.3.8 and 

2.4.9. 

Axial Turbine Inputs-Ammoni a 

if Stages it e act ion RPM Units Feasible 

[ 1 m 

2 .5 900 60 4 Yes 
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# Stages 



Reaction 



RPM C<, 
lm 



Units Feasible 



10 


.5 - 


720 


60 


4 


No 


5 


.5 


720 


60 


4 


Yes 


2 


.5 


7-20 


60 


4 


Yes 


10 


.5 


900 


60 


2 


No 


5 


.5 


900 


60 


2 


No 


2 


.5 


900 


60 


2 


Yes 


10 


.5 


720 


60 


2 


No 


5 


.5 


720 


60 


2 


No 


2 


.5 


720 


60 


2 


Yes 


5 


.5 


900 


60 


4 


Yes 


10 


.5 


900 


60 


4 


No 


2 


.5 


720 


60 


10 


Yes 


5 


.5 


720 


60 


10 


Yes 


10 


.5 


720 


60 


10 


Yes 


2 


.5 


900 


60 


10 


Yes 


5 


.5 


900 


60 


10 


Yes 


10 


.5 


900 


60 


10 


No 


5 


.5 


900 


60 


4 


Yes 


5 


. 6 


900 


60 


4 


Yes 


5 


.7 


900 


60 


4 


Yes 


5 


.8 


900 


60 


4 


Yes 


5 


.5 


900 


70 


4 


No 


5 


.6 


900 


70 


4 


No 


5 


.7 


900 


70 


4 


Yes 


5 


.8 


900 


70 


4 


Yes 


5 


.5 


900 


80 


4 


No 


5 


. 6 


900 


80 


4 


No 


5 


.7 


900 


80 


4 


No 


5 


.8 


900 


80 


4 


No 


20 


.2 


720 


60 


20 


No 


20 


.2 


600 


60 


20 


No 


10 


.2 


600 


60 


20 


Yes 


20 


.2 


600 


50 


20 


No 


10 


.3 


600 


60 


20 


Yes 


20 


.3 


600 


60 


20 


No 


20 


.3 


600 


50 


20 


Yes 


20 


.3 


600 


40 


20 


Yes 


20 


.3 


600 


30 


20 


Yes 


20. 


.3 


600 


20 


20 


Yes 




Axial 


Turbine 


Inputs-Pr eon 


21 




2 


.7 


600 


60 


20 


Yes 


5 


.7 


600 


60 


20 


No 


10 


.7 


600 


60 


20 


No 


2 


.7 


720 


60 


20 


No 


5 


.7 


720 


60 


20 


No 
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# Stages 


Reaction 


RPM 


e< 

lm 


Units 


Feasible 


10 


.7 


720 


60 


20 


No 


2 


.7 


720 


70 


20 


No 


5 


.7 


720 


70 


20 


No 


10 


.7 


720 


70 


20 


No 


2 


.7 


720 


80 


20 


No 


5 


.7 


720 


80 


20 


No 


10 


.7 


720 


80 


20 


No 


2 


.2 


720 


60 


20 


No 


2 


.3 


720 


. 60 


20 


No 


2 


.4 


720 


60 


20 


No 


2 


.6 


720 


60 


20 


No 


2 


.7 


720 


60 


20 


No 


2 


.8 


720 


60 


20 


Yes 


2 


.9 


720 


60 


20 


No 


2 


.7 


600 


70 


20 


No 


5 


.7 


600 


70 


20 


No 


10 


.7 


600 


70 


20 


No 


2 


.7 


600 


80 


20 


No 


5 


.7 


600 


80 


20 


No 


10 


.7 


600 


80 


20 


No 




Radial 


Turbine 


Input s-Ammonia 




RPM 


Units 




# Blades 




Feasible 


3600 


13 




12 




Yes 


3600 


12 




12 




Yes 


3600 


11 




12 




Yes 


3600 


10 




12 




Yes 


3600 


9 




12 




Yes 


3600 


8 




12 




Yes 


3600 


7 




12 




Yes 


3600 


6 




12 




Yes 


3600 


5 




12 




Yes 


3600 


4 




12 




No 


3600 


14 




12 




Yes 


3600 


15 




12 




Yes 


3600 


16 




12 




Yes 


3600 ' 


17 




12 




Yes 


3600 


18 




12 




Yes 


3600 


19 




12 




Yes 


3600 


20 




12 




Yes 


3600 


21 




12 




Yes 


3600 


22 




12 




Yes 


3600 


23 




12 




Yes 


3600 


24 




12 




Yes 


1800 


1 




12 




No 
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RPM 




Units 


# 


Blades 


Feasible 


1800 




2 




12 


Yes 


1800 




3 




12 


Yes 


1800 




4 




12 


Yes 


1800 




5 




12 


Yes 


1800 




6 




12 


Yes 


1200 




1 




12 


Yes 


1200 




2 




12 


Yes 


900 




1 




12 


Yes 






Uadial Turbine 


Input s-Fr eon 21 




3600 




250 




15 


No 


3600 




500 




15 


Yes 


3600 




750 




15 


Yes 


3600 




1000 




15 


Yes 


3600 




1250 




15 


Yes 


3600 




1500 




15 


Yes 


3600 




1750 




15 


Yes 


3600 




2000 




15 


Yes 


3600 




2250 




15 


Yes 


3600 




2500 




15 


Yes 


3600 




2750 




15 


Yes 


3600 




3000 




15 


Yes 


3600 




3250 




15 


Yes 


3600 




3500 




15 


Yes 


3600 




3750 




15 


Yes 


3600 




4000 




15 


Yes 


3600 




4250 




15 


Yes 


3600 




4500 




15 


Yes 


3600 




4750 




15 


Yes 


3600 




5000 




15 


Yes 


3600 




5250 




15 


Yes 


3600 




5500 




15 


Yes 


3600 




5750 




15 


Yes 


3600 




6000 




15 


Yes 


3600 




6250 




15 


Yes 


P.4 


Sample 


Output 










The o u 


tput for both working fluids produced 


many designs 


which 


could 


be called feasible. 


A 


representative 


sample of 



axial and radial designs for each fluid is included on the 
following pages. These designs are not necessarily optimums, 
as the object of the runs was to test the program rather 
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than optimize a design. Data from these runs was used to 
substantiate the evaluation of the model, sections 3.5 and 
3.6. These examples do, hoVever, illustrate several points 
about the Zener Cycle and will be the basis for the next 
section . 

Axial Turbine-Ammonia 

Assumed 5. Stages 0.50 Reaction at mean 720. RPM 
Assumed 60.00 ALFAMN 4. Units 6050.00 m 

Last stage tip rei Mach no. 0.306 Last stage hub real Mach 

no . 0.314 

First stage tip diam. 8.51 Blade height 0.9599 if blades 50. 

Eff. 0.9589 

Last stage tip diam. 8.84 Blade height 1.2949 ft blades 38. 

Eff. 0.9598 

Machine efficiency 0.9594 
Max cent, stress 3003.50 
Unit length 2.492 

Aspect ratio, first stage 2.17 Last stage 2.34 

Velocity Triangles 

First Stage 

Hub 



UHN 


248.3 


VTHN 325.9 


VX 164.2 VHN 


365 .0 


VHN 181.7 


WHR 


297.7 


ALFAHN 63.3 


BETAliN 25.3 


BETAHR 


-56.5 


Mean 










UMN 


284.5 


VTMN 284.5 


VX 164.2 V MN 


328.5 


VMN 164.2 


WMR 


328.5 


ALFAMN 60.0 


BETAMN 0.0 BET AMR - 


•60.0 


Tip 












UTN 


320.6 


VTTN 252.4 


VX 164.2 VTN 


301.1 


VTN 177.9 


WTR 


360.3 


ALFATN 56.9 


BETATN -22.6 


BETATR 


: -62.9 


Last 


. Stage 










Hub 












UiiN 


235 .6 


VTHN 343.4 


VX 164.2 VHN 


380.6 


VHN 196.4 


WHR 


287.2 


ALFAHN 64.4 


BETAHN 33.3 


BETAHR 


-55.1 


Mean 










UMN 


284.5 


VTMN 284.5 


VX 164.2 VMN 


328.5 


VMN 164.2 


WMR 


328.5 


ALFAMN 60.0 


BETAMN 0.0 LET AMR - 


60.0 


1 X ]) 

UTN 


333.3 


VTTN 242.8 


VX 164.2 VTN 


293.1 


VTN 187.5 


WTR 


371.5 


ALFATN 55.9 


BETATN -28.9 


BETATli 


. -63.8 
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Axial Turbine-Freon 21 



Assumed 2. Stages 0.70 Reaction at mean 720. RPM 
Assumed 60.00 ALFAMN 20. Units 30220.00 m 



Last stage tip real Mach no. 0.741 Last stage hub real 

Mach no. 0.542 

First stage tip diam. 8.87 Blade height 1.9793 ji blades 84 

Eff. 0.9712 

Last stage tip diam. 9.67 Blade height 2.7744 jf blades 23. 

Eff. 0.9701 

Machine efficiency 0.9706 
Max cent, stress 5878.27 
Unit length 0.802 

Aspect ratio, first stage 10.00 Last stage 4.59 



Velocity Triangles 



First Stag 
Hub 

UHN 185.2 
WHR 205.9 
Mean 

UMN 259.8 
WMU 275.0 
Tip 

UTN 334.4 
WTK 346.3 
Last Stage 
Hub 

UHN 155.2 
WHR 179.4 
Mean 

UMN 259.8 
WMR 275.0 
Tip 

UTN 364.4 
WTU 375.4 



VTHN 218.7 
ALFA1IN 67.6 

VTMN 155.9 
ALFAMN 60.0 

VTTN 121.1 
ALFATN 53.4 



VTHN 260.9 
ALFAHN 71.0 

VTMN 155.9 
ALFAMN 60.0 



VX 90.0 VHN 236.5 WHN 96.0 
BETAI1N 20.4 BETAHR -64.1 

VX 90.0 VHN 180.0 VsMN 137.5 
BETAMN -49.1 BET AMR -70.9 

VX 90.0 VTN 150.9 VTN 231.5 
BETATN -67.1 BETATR -74.9 



VX 90.0 VHN 276.0 WHN 138.8 
BETAHN 49.6 BETAHR -59.9 

VX 90.0 VMN 180.0 WHN 137.5 
BETAMN -49.1 BET AMR -70.9 



VTTN 111.2 VX 90.0 VTN 143.0 WTN 268.8 
ALFATN 51.0 BETATN -70.4 BETATR -76.1 



Radial inflow Turbine-Ammonia 



Assumed 3600 RPM 17. Units 12. Blades 6050.00 WW 

Nozzle inlet diam. 6.99 Blade height 0.2338 
Nozzle outlet diam. 4.63 

Rotor inlet diam. 3.70 Outlet tip diam. 2.91 

Rotor hub diam. 1.52 

Unit length 1.09 

Max scroll area 3.972 

Specific speed 0.127727 

Efficiency 0.9398 

Velocity Triangles 
Nozzle inlet 

VI 362.8 VT1 229.8 VM1 280.7 ALFA1 39.3 
Rotor Inlet 

U3 696.8 VT3 580.7 VM3 364.9 V3 685.8 W3 382.9 
ALEA3 57.9 BETA3 -17.7 
Rotor Outlet 
Hub 

U4II 287.4 V4H 241.4 W4H 375.3 BETA4II -50.0 
Mean 

U4M 418.1 V4M 241.4 W4M 482.8 BETA4M -60.0 
Tip 

U4T 548.8 V4T 241.4 W4T 599.5 BETA4T -66.3 



Radial Inflow Turbine-Freon 21 



Assumed 3600 RPM 1250. Units 15. Blades 30220.00 WW 



Nozzle inlet diarn. 2.75 Blade height 0.1046 
Nozzle outlet diain. 2.04 

Rotor inlet diam. 1.62 outlet tip diam. 1.23 

Rotor hub diam. 0.72 

Unit length 0.45 

Max scroll area 0.628 

Specific speed 0.113886 

Efficiency 0.9345 



Velocity Triangles 
Nozzle Inlet 

VI 173.4 VT1 124.8 VM1 120.4 
Rotor Inlet 

U3 305.8 VT3 265.0 VM3 128.1 
ALEA3 64.2 BETA3 -17.7 
Rotor Outlet 
Hub 

U4H 135.1 V4H 105.9 W4H 171.7 
Mean 

U4M 183.5 V4M 105.9 W4M 211.8 
Tip 

U4T 231.8 V4T 105.9 V4T 254.9 



ALFA1 46.0 
V3 294.3 W3 134.4 

BETA4H -51.9 
BETA4M -60.0 
BETA4T -65.4 
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F.5 Comment on the Zener Cycle 

Using the sample output shown in the previous sections, 
it is possible to calculate a volume required for each machine. 
Multiplying this by the number of units gives a volume to 

produce 100,000 kilowatts of power. This approximate volume 

3 3 

is 600 ft for Axial with Ammonia, 1200 fv for Axial with 

Freon 21, 700 ft^ for Radial with Ammonia, 3700 ft^ for 

Radial with Freon 21. This is interesting in that it shows 

ammonia to be the more advantageous of the two fluids in 

terms of size. This would be expected, however, in that for 

the cycle points defined, ammonia has a greater enthalpy 

change per pound than Freon. Freon has the advantage of 

being safer to work with, however. 

A second fact is apparent in the volume consideration 
in that the estimated volume was in the order of one to two 
thousand cubic feet. Zener estimated a turbine module to be 
8 x 8 x' 40 feet or about 2560 cubic feet. This shows that 
Zener’s estimation was reasonable for plants of 100,000 
kilowatt output. Assuming this size plant, then to supply 
the output of 60 billion kilowatts mentioned in section F.l, 
600,000 individual plants would be required. This also 
assumes that the heat exchangers are of comparable size. 

(Jompar i sons could be made for the amount of working 
medium required, amount of cooling water and others. That 
is not really the point of the evaluation, however. Exami- 



nation of the sample designs does show, for at least the 
turbines, that the proposal is at least within the realm 
possibility if not reasonability. 
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